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Abstract 

The quantum mechanical formula for Mayer's second cluster integral for the 
gas of relativistic particles with hard-core interaction is derived. The proper 
pion volume calculated with quantum mechanical formula is found to be an 
order of magnitude larger than its classical evaluation. 

The second cluster integral for the pion gas is calculated in quantum me- 
chanical approach with account for both attractive and hard-core repulsive 
interactions. It is shown that, in the second cluster approximation, the repul- 
sive 7T7r-interactions as well as the finite width of resonances give important 
but almost canceling contributions. In contrast, an appreciable deviation from 
the ideal gas of pions and pion resonances is observed beyond the second clus- 
ter approximation in the framework of the Van der Waals excluded-volume 
model. 

I. INTRODUCTION 

Thermal models ("fireball") have been popular for decades (see, e.g. |L|) to fit the data on 
multiparticle production in high energy nucleus-nucleus collisions (see, e.g. [@ and references 
therein). The ideal gas (IG) model of noninteracting hadrons and resonances which has 
mostly been empoyed to extract temperature T, baryonic chemical potential hb, etc., from 
fits to the data is however not adequate for this purpose. This is among other things because 
of the following two reasons: 

• The ideal gas model ignores the finite width of the resonances, while most of them 
have a width comparable to or even larger than the typical temperatures of the hadron 
gas 120 -T- 180 MeV. This leads to underestimation of the attraction between hadrons. 

• The IG model does not take into account nonresonance interaction between hadrons, 
in particular the repulsion. As a result, in the description of the hadron yields data 
for the AGS and SPS energies the IG model leads to artificially large particle number 
densities, e.g. p = 1.25 fm" 3 at T = 185 MeV and fi B = 270 MeV [|], which hardly 
can be consistent with a picture of a gas of point-like, noninteracting hadrons. 

The procedure of Ref. [|J introduces the Breit-Wigner mass spectrum of resonances. 
However, this widely used procedure works for narrow resonances only. It was found that 
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it is insufficient in the realistic case || as it does not take into account correlated particle 
pairs appearing along with resonances in the hadron gas. Therefore, the standard procedure 
underestimates tha attractive part of hadron interactions. 

To solve the second problem the procedure, which allows to take into acount finite 
particcle volume, was proposed by Hagedorn and Rafelski ||. The excluded- volume Van der 
Waals equation of state was derived in Ref. [0] and used by several authors (see, e.g. @,|3||| 
and references therein). Recently, this procedure was generalized to multicomponent M and 
relativistic particle systems Still, the proper particle volume was so far calculated by 
classical statistical mechanics formulae. 

Our aim is to calculate Mayer's cluster integrals (CIs) for the hadron gas from the 
available data on the hadron scatterings using correct quantum formulae and use them for 
fixing the parameters of the Van der Waals excluded volume model. In the present paper a 
first step in this direction is made, namely we calculate the second cluster integral (CI) in 
the case of a pure pion gas for a wide temperature range and consider the contribution of 
the repulsive part of the nir interactions into the CI as an excluded volume of the Van der 
Waals model. 

The article is organized as follows: in section II we derive the formula for the second 
cluster integral taking into account relativistic effects as well as the isospin of the pion. 
Section III is devoted to the hard-core repulsion at the quantum level. The domain of 
applicability of the classical formulae is found. The resonance attraction will be considered 
in section IV. The conditions which allow to use the narrow resonance approximation (NRA) 
and the Bright-Wigner formula of Ref. |4| will be studied. In section V the CI for the 
pion gas is calculated from the experimental data on the 7T7r-scattering. The results are 
compared with various approximations. In section VI the interacting pion gas is studied 
in the framework of the excluded-volume Van der Waals approach. The conclusions and a 
discussion of the results are given in section VII. 



II. GENERAL FORMULAE 

Quantum mechanical formula for a calculation of Mayer's second cluster integral 6 2 in 
the case of nonrelativistic zero isospin I = particles was considered in Ref. JTlJ (see also 
T2J). The pions, however, have nonzero isospin I = 1 and the temperature of interest can 



be comparable to or larger than the pion mass. Therefore, for an adequate description of 
the pion gas a generalization of the formulae given in Refs. [ |TT| , |T2] for relativistic particles 
carrying nonzero isospin is needed. 

We start from the canonical partition function for N identical particles in the volume V 
at the temperature T 

Z{V, T,N) = J d 3N r K(*i, r 2, • • • r N ) exp (-^) * a (n, r 2 , . . . r*), (1) 

where H is the Hamiltonian operator and {^ a } is a complete set of orthonormal wave 
functions in co-ordinate representation. With the notations 

W n (t u t 2 ,...t n ) = N\ J2 K fa, r 2 , ■ • • *n) exp ( -— ) tt a (n, r 2 , . . . r N ) (2) 

a 



T 
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one gets 

Z(V, T, N) = ^ J d 3N rW N ( ri , r 2 , . . . r N ) . (3) 
The function Wi(ri) can be calculated in the thermodynamical limit V — > oo 

^ (ri )^_exp(--) (4, 



p,tj 



< 2/ » +1 >7($F exp (- 



Vp 2 + m 2 
T 



g<p(T;m) 



where J and m are, respectively, the particle isospin and mass (tj = —I ,...,+I is the 
isospin projection), g = (2I + 1) is the isospin degeneration factor[| and <f>(T; m) can be 
expressed via K 2 modified Bessel function 

1 f°° 2 / Vp 2 + m 2 \ m 2 T (m\ 
mm) = — 2 J Q pdpex^ —j = -^K^-) . (5) 

The asymptotics of 4>(T; m) in the nonrelativistic, m >> T, and ultra-relativistic, m << T, 
limits are 

f i^f 2 exp(-m/T) , m»T 
ct>{T-m) ~ (6) 

Following Refs. JT^JT^] we introduce the functions C//(r 1; r 2 , . . . r^): 

W2( ri ,r 2 ) = C/^rOC/i^) + C/ 2 (ri,r 2 ) (7) 
^ 3 (ri,r 2 ,r 3 ) = ^i(r 1 )f/ 1 (r 2 )[/ 1 (r3) + f/i(r 1 )f/ 2 (r 2 ,r 3 ) 

+^i(r 2 )C/ 2 (r 3 , r x ) + l/ifo^fa, r 2 ) + E/gfa, r 2 , r 3 ) 

etc. 

and define Mayer's CIs^ 



If not only the isospin but also the spin has a nonzero value, Jo, the degeneration factor has the 
form g = (2J + 1)(2J + 1). 

2 The normalizations of the CIs in Eq. (|8|) are different from that of Ref. |]l^] and correspond to 
the definition used in [|rj[]. The CIs (||) have dimensionality [volume]' -1 , while in Ref. (T^] CIs are 
dimensionless. 
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It is easy to see that 



h = 1 . 



(9) 



Substituting the expression (|7|) for the function W 2 into Eq. ([3]) one gets the two-particle 
partition function expressed via the CIs 



Z(V,T,2) = g 2 



MV, T)0(T; m)Vf + b 2 (V, T)0 2 (T; m)V 



(10) 



For arbitrary iV the expression of the partition function via bi reads [12 



N 



Z(V,T,N) = Y.W—AHV^ig^m^Vr 
{m l} i=i ™< ! 



(11) 



where the sum runs over all sets of nonnegative integer numbers {mi} satisfying the condition 



N 



Y, lm i = N ■ 



(12) 



i=i 



Introducing the absolute activity 0], which in the relativistic case takes the form 

z = 9<f>{T]m)exv(Jpj , (13) 
where /i is the chemical potential, the grand canonical partition function has the form 



Z(V,T,fjL) = ^2exp(^-)Z(V,T,N) = exp(V$>(F,T); 



(14) 

N=l V -i / \ 1=1 / 

From Eq. ([14]) one can find the cluster expansion of the pressure and the particle density: 

log Z(V,T,fi) 



p(T,li) = T \im 

V — >oo 



V 



n(T, fx) 



. T d log Z(V,T, li) 



i=i 

oo 



lim 



dfx 



where 



k{T) = lim bi{V,T). 



V^oo 



(15) 
(16) 

(17) 



Substituting the particle density (|16"D into the virial expansion^ for the pressure 



3 Sometimes the term 'virial expansion' is used in place of 'cluster expansion' |H,22]. We prefer 
to use the standard terminology [[1^,14.15]: 'cluster expansion' for the expansion in powers of the 



activity and 'virial expansion' for that in powers of the particle density. 
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p(T, n) = Tj2 



(18) 



i=l 



and equating the coefficient of each power of z with Eq.(|15]) one obtains the following 
expressions for the virial coefficients [12j in terms of the CIs: 



a% = 1 

a 2 = —b 2 

a 3 = Ab\ - 2b 3 

a 4 = -20b 3 2 + 18& 2 &3 - 36 4 



(19) 



Let us represent the CIs as a sum of two terms: 



where &j are the CIs for the IG and b\ l> appear due to the particle interaction. In the classical 



bf + b? 

■,(0 



I > 1, 



(20) 



(Boltzmann) gas one obtains b\ = for all I > 1. In the quantum case, b\ u> are nonzero due 
to Bose (Fermi) effects and can be easily found for arbitrary /. For noninteracting particles 
the logarithm of the expression ( JPty should coincide with the well-known expression for the 
logarithm of the ideal gas grand canonical partition function 



l(0) 



log Z^(V,T,fi) = ± 9 V 



d 3 p 
(2^)3 



log 



1 ± exp 



yU — y/p' 2 + m' 2 

f 



(21) 



(the upper (lower) sign corresponds to Fermi-Dirac (Bose-Einstein) statistics). One can 
expand the logarithm in the integrand and perform the integration 



log Z^(V,T,fi) = 9 VJ2 



(Tl) 



i+i 



i=i 

oo 



d 3 p 



exp 



I (fi — \?p 2 + rn? 



gV }^ — J — ex P y-f J <f>{T/l; m) . 



Comparing the last expression with Eq. ([14]) gives 



(22) 



6f> 



(Tl) i+1 0(T//;m) 



In the nonrelativistic limit Eq.( 



Ig*- 1 [<KT;m)y 
is reduced to 



(23) 



6[°) = { T i) l+ H^' 2 
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(24) 



where A is the thermal wave length 



A = 



(25) 



The expression (|24j ) coincides for Iq = with the corresponding formulae of Ref. |12| (up to 
the dimensional factor A 3 *-'^, because of different normalization in Eq. @). 

Using Eqs. (|I0| ) and (|20|) one can express b$ via differences of the two-particle partition 
functions for real and ideal gases: 

( „ Z(\/,T,2)-ZW(\/,T,2) 

2 " V^(T;m)]« W 

Let us calculate Z(V,T, 2). A complete set of the orthonormal state vectors in the two 
particle system can be constructed from the following wave functions 

e i(P,R) 

| a) = |P,a) = — 7^ | a) , (27) 



V 

where P is the total momentum of the system, R is the radius-vector of its center of mass, 
and | a) form a complete set of orthonormal state vectors of the system in the center of mass 
frame (c.m.f.) satisfying the Schrodinger equation 

H \&) = £& 1 5) (28) 

with the normalization condition 

(a'\a) = 5 & 'a • (29) 
The wave function ([27D thus satisfies the following equations 



//|P,a) = v /P 2 + 4|P,a), (30) 

(P',a'|P; a ) = <Wp'p • (31) 
The expression for Z(V, T, 2) in terms of the introduced wave functions has the form 



Z(V,T,2) = £( P ,a|ex P (-f)|P,a) = E ex p(-^^| • (32) 

P,5 \ J- / P,a \ 1 / 

In the thermodynamical limit V — > oo the summation over P can be replaced by the 
integration and one finds 



, Mp I /p2 i £ 2 \ 

Z(V,T,2) = ^Vj ^exp 1- V — a \ = V^me.) ■ (33) 

The states of two spinless particles in their c.m.f. can be enumerated by the following 
quantum numbers: the radial momentum q (or, alternatively, the energy in c.m.f. e(q) = 
\J q 1 + m 2 ), the orbital angular momentum L, its projection tul, the total isospin / and its 
projection ti, e.g. a = (q, L,rriL, I,ti). Eq. (|33|) can be rewritten explicitly 



l We assume that the particles do not form bound states with energy e < 2m. 
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2I I L -.qc 

Z(V,T,2) = E E E / dqg LmLltl (q)me(q)) , (34) 

7=0*1 = -/ L m L =-L J ° 

where gtm L it I (q) is the density of states with the given set of quantum numbers. The sum 
J2' extends over only those values of L that satisfy the symmetry properties of the wave 
function. For spinless bosons it takes even value if the isospin part of the wave function is 
symmetric and odd values if it is antisymmetric. In the case of integer isospin particles like 
pions this means Q 

J 0,2,4, 6, ... for even I , . 

\ 1,3,5,7,... forodd/ ' [6b) 

Substituting Eq.fl3"4"|) into (|2T)| ) one gets the following expression for the second CI 

$ = T-T^pE E E' E / dq(g LmLlt M-9 ( l Llt Mme( q )), (36) 

where gLm L itz * s ^ ne state density for the IG. The difference ghmLiuW) — 9i}n L itX c i) i n the 
thermodynamical limit can be expressed via phase shifts of two-particle scattering 5L mL i tl (q) 

OB ' 

/ \ (o) / x 1 d5 Lm It {q) 

9Lm L ItM) ~ 9Lrn L ItM) = ~ ^ • ( 37 ) 

Using the expression for <f>(T] m) one gets 



&w = - V y y Y dqe\q) LmL " iyH, K 2 {e{q)/T) . (38) 

In the case of hadron gas the phase shift does not depend on the angular momentum pro- 
jection mi (no external fields) and the isospin projection tj (if only strong interactions are 
taken into account). This simplifies the last formula: 



(39) 



Performing a partial integration and taking into account the properties of the Bessel func- 
tions, one gets 

4" = m *T> [9 K 2(m/ T W § P 2 ' + W + £ ^''^ (?) ' <40) 
In the nonrelativistic limit the formula (|3~9T) is reduced to 

,2 



which again at Iq — coincides with the corresponding formulae of Refs. jlllJEql , up to the 
factor A 3 . 



5 In the case of half-odd isospin scalar bosons we would have the opposite rule: even L for odd I 
and odd L for even /. 
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III. HARD CORE REPULSION 



The hard core repulsion plays an important role in the phenomenological description 
of the 7T7r-scattering: the phase shift data for the isospin state 1 = 2 can be successfully 
described assuming hard core repulsion between two particles |L6| . The best fit of the phase 
shift in So-state (I = L = 0) can be obtained by assuming hard core repulsion in addition 



to resonance attraction |T^,|T8[. Hence we start our analysis from applying Eqs.(|3^-^1|) to 
hard sphere model.Q 

The radial part of the wave function in the c.m.f. of two particles interacting by hard core 
potential in the state with orbital momentum L and radial momentum q can be represented 
in the following way 

, , v _ ( for r < r , , 

1 C(cos5l + sin<5k yi(gr)) for r > r ' 

where r is the distance between particle centers, r$ is the minimal admitted value for r 
(that is the doubled radius of the particle considered as a hard sphere), j n (z) and y n (z) are 
spherical Bessel functions, C is the normalization constant and 8l is fixed by the condition 

cos 8 L j L {qr ) + sin 8 L y L (qr ) = . (43) 

Using asymptotic properties of the spherical Bessel functions, it is easy to see that 8l has a 
meaning of the phase shift describing scattering of two hard spheres: 



sin (qr - l -f + 8 A 



qr 

The derivative of the phase shift is found to be 



oo . (44) 



dh = ± arctan ( 3L{qrp) \ = tq ^ 

dq dq \yL,(qr ) J (qr ) 2 [jUqr ) + yl(<Fo)] ' 

where the formula for the Wronskian [19| 

W(j L (z),y L (z)) = z- 2 (46) 

has been used. 

Let us introduce the functions 

K±{Z)= ^z*\jl(z) + yl(z)] ' (47) 

where the sum J2l runs over either even (superscript '+') or odd (superscript ' — ') nonneg- 
ative numbers. 



6 Relativistic consideration of a hard sphere model by no means can be consistent. Still, as far as 
this model describes experimental data on low energy 7T7r-scattering, we find it phenomenologically 
satisfactory. 



S 



Expanding k around zero, one gets 



K + {z) ~ 1 + + 0{Z 6 ) , 
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K -( z )^3 z *-3z 4 + ^z 6 + 0(z s 



(48) 
(49) 



It has been checked numerically that the asymptotic behavior of k (z) at large z with a 
high accuracy can be presented by the formula 



« ± (z) ~ ^ 2 + ^z + i + 0(2- 1 ) . 



(50) 



The CI ( P^D for the case of hard sphere model can be represented in the formfj 



6 (0 

°2 



1 2/ 

-^(2/ + l)6 2 ± (r 0) T, 



m 



(51) 



where superscript '+'(' — ') corresponds to even (odd) values^ of / and bf(r ) is expressed 
via the function k (z) 



bf 



lure 



m*T[K 2 (m/T)} 2 Jo 



dq[e(q)} 2 K ± (qr )K 2 



T 



(52) 



In the nonrelativistic approximation the expression for bf(r ) is reduced to 

2y / 2 r°° ' ' ' 

bf(r ,T,m) = A 3 ro / dqK ± (qr ) exp 

7T JO 



mT 



(53) 



Substituting Eq. (E8j) into the last formula one gets the nonrelativistic expression for 



b 2 (r ,T, m) at small r (r << A) 



6+(r , T, m) ~ -2A 2 r (l + -7r 2 (r /A) 4 + O [(r /A) 6 
62 (r„, T, m) ~ -6vrr 3 (l - 37r(r /A) 2 + ^7r 2 (r /A) 4 + O [(r /A) 



(54) 
(55) 



in the opposite case r >> A the nonrelativistic expression for 6^(r ,T, m) can be obtained 
using Eq. 



bf(r ,T,m) = --vrr 3 , 1 + 



3^2 A j^A 2 
4 r 27r rg 



+ 



A 

Jo. 



(56) 



In the simplest version of the hard sphere model when the core radius tq does not depend on 
/, the formula (]5l|) could be reduced to = bf(ro,T,m). However, in the case of realistic -kit 
interaction every isospin state has its own Tq{I) |fl7| , ^6| , ^8| . 

8 Again, for the case of half-odd isospin particles the opposite rule would be valid. 
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The pion thermal wave length at the temperature T = 50 -j- 200 MeV ranges between 



3-^6 fm. From Eq. ( pq) one sees that the classical formula [15 



6f(r ,T,m)^-^7rr 3 (57) 

(the particle volume multiplied by 4) would give a reasonable approximation only at unreal- 
istically large hard core radius ro > 50 fm. The hard core radii found from the 7T7r-scattering 



are much smaller: ro = 0.60 fm in the Sq state ]IS| and ro = 0.17 fm at I = 2 [[UJ. (No 
evidence of hard core repulsion was found in Pi-state (I=L=1)). In this case the value of 
6^"(ro,T, m) can be estimated from formula (|54|) . The results are presented in Fig.II]. It is 
seen that in contrast to the classical case the quantum treatment leads to a rather strong 
dependence of the second CI on the temperature (approximately proportional to 1/T) even 
in the nonrelativistic approximation. Numerical calculations show that relativistic effects 
make this dependence even stronger and essentially reduce at high temperature the CI with 
respect to its nonrelativistic value (see Fig. |l]). Both relativistic and nonrelativistic quan- 
tum formulae give a much (1-^2 orders of magnitude) larger value than those given by the 
classical formula 0.45 fm 3 and 0.01 fm 3 for r = 0.60 fm r = 0.17 fm, respectively. 

As can be seen from Fig.[l], the relativistic effects cannot be ignored even at relatively 
low temperatures. Therefore, only the relativistic formula ([52]) is used in the following 
calculations. 



IV. RESONANCE ATTRACTION 

The phase shifts of hit elastic scattering can be approximately described by assuming 
that the attractive parts of the interaction appear due to the propagation of resonances in 
the s-channel of the reaction (see [16|- |T8| , p0f and references therein). 



The distinctive feature of the resonance interaction is the rapid growth of the phase shift 
by it radian in the vicinity of the pion momentum q r , which is related to the resonance mass 



M r = 2^ q 2 r + m? . (58) 

In the limit r r — > (T r is the resonance width) the derivative of the phase shift can be 
approximated by the Dirac delta-function: 

— — — «7r%-g r ) . (59) 

In this approximation, which we will call the 'narrow resonance approximation' (NRA), the 
expression (|39"D for b% yields: 



2 ^(2/ r + l)(2L r + l) ^ ' r 2 , (60) 
9r mT-m)Y 

where I r , L r and M r are, respectively, the resonance's isospin, spin and mass, with the index 
r running over all resonances in the two-pion system. 
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Eg. fl60D allows to rewrite the expression for the grand canonical partition function (|i"4"|) 
in the following form 



Z(V,T,/j) = exp 



\l=l r 



where 



:,. = <>{T: M r ) ox]> ( ^ 



(61) 



(62) 



is the absolute activity of the resonance r with degeneration factor g r = {21 r + l)(2L r + 1). 
The expression (|61| ) is nothing else than the partition function for a mixture of ideal gases 
of pions and two-pion resonances^. This recovers the well known result of Ref. j^TJ that 
narrow resonances contribute to the partition function as an ideal gas of stable particles. 



The quantitative criterion for an applicability of the NRA was found to be fl2T 

T r «T. 



(63) 



The resonances appearing in 7T7r-scattering do not satisfy this criterion: most of them (p(770), 
/o(980), /2(1270) , p3(1690)) have widths comparable with a typical temperature of the 
hadron gas and the width of /o(400 — 1200) (known also as the a) is a few times larger 
then the temperature. Therefore, it is necessary to take into account the finite width of the 
resonances. 

The scalar- isoscalar resonances contribution to the 7T7r-phase shift in the Sq state can be 
parametrized in the following way |I7 |: 



q M 2 

tan5 r (g) - 



q T e{q)M 2 -e 2 {q) 



r = a, / (980) 



(64) 



For parametrization of nonzero (iso-)spin resonances we shall use the following formula [20 



tan5 r (g) 



' q\ r+ M r x r T r D Lr (q r R r 
^qr) M?-e 2 (q) D Lr (qR r ) 



p(770), / 2 (1270), p 3 (1690), (65) 



where x r is the inelasticity, i.e. the decay fraction of the resonance into two pions, R r is the 
so-called interaction radius and the functions Dl(z) are given by the formulae 



D 1 (z) 
Ds(z) 



1 + z 2 

9 + 3z 2 + z A 

225 + 45z 2 + 6z 4 + z 6 



(66) 



The resonance parameters are given in the Table |. 



9 The fact that the resonance gases are classical is an artifact of the second cluster approximation. 
The quantum correction would appear from the 4-th and higher cluster integrals for the interacting 
pions. 
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It is easy to see that if a resonance lies far from the threshold: 



M r - 2m » T r (67) 

both formulae are reduced to 

- wh < 68 > 

In this case the activity of the resonance can be represented in the form 

2I 



E* = EE' [°° d <( £ )( 21 + !)( 2L + WW £ ) ex p {%) > 

where the resonance profile function is given by the Breit-Wigner formula: 



(69) 



i r r 

C(£) = 2^(e- M r ) 2 + (r r /2) 2 (70) 

From this we conclude that the procedure of Ref. |IJ where the profile function was postulated 
to be 

C(e) = — (71) 

U ; (e - M r ) 2 + (r r /2) 2 1 ; 

with normalization constant £ fixed by the condition 

2m & (e-M r ) 2 + (r r /2) 2 = 1 (72) 



becomes valid in the limit fl67Q . The a-resonance obviously does not satisfy this condition, 
even for the p(770) the difference M r — 2m is only about 3 times larger then the width. We 
have calculated the contributions of these two resonances into the second cluster integral 
&2 using the parametrizations (0) and ( |65D and compare them with the corresponding 
approximate values found in the framework of the procedure of Ref. M : 

^ = m*T[gK 2 (m/T)}^ I dq£2{q) (a - M r fl {T r /2f ^T* 2 (t) ' (?3) 

and in the NRA (|60|) . The results are shown in Figs. ^| and [3]. As can be seen from Fig.|2|, 
for a-resonance the both approximations essentially underestimate the CI. It is interesting 
to mention that at T > 150 MeV the formula of Ref. 0] gives slightly worse result than even 
that of the NRA. In the case of the p-resonance this formula systematically overestimates 
the CI in contrast to the a case. The role of the resonance width becomes small at high 
temperatures and all three formulae give comparable results in both (a and p) cases. 
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V. INTERACTING PION GAS 



Both type of interaction: hard core repulsion and resonance attraction are present in 
the pion gas. The phase shift for 7T7r-scattering in the So state at the center of mass energy 
below 1 GeV can be represented as a sum of three terms |T^^8|: 



M<?) = $a(q) + S fo (q) + 6 BG (q) . (74) 

The background term Sbg is related to the hard core repulsion 

W<?) = ~r q , r = 3.03 GeV" 1 (75) 

and two first terms describe attraction due to the resonances a and /o an d are parametrized 
by the formula (|B"4]). The contributions of the Sq state to the CI are shown in FigfJ. The 
attractive part is larger in absolute value than the repulsive part, so that total contribution 
of the S state is positive. 

The interaction in the 5*2 state has a purely repulsive nature (no exotic resonances with 
isospin 1 = 2 have been found). The phase shift can be successfully fitted by the hard core 
formula: 

M?) = ~r q , r = 0.87 GeV" 1 . (76) 

As it is seen from FigfJ the absolute value of the negative contribution of S2 state into the 
CI is slightly larger than that of the positive contribution of the Sq state, so that these two 
quantities almost cancel each other. The resulting contribution of the S-state into CI is 
negative and a few times smaller than those of the So- and S2-states separately. Due to the 
small value of the total S-state contribution into CI the Pi state becomes important already 
at relatively low temperature. 

The phase shifts of 7T7r-scattering in Pi, D2 and Pi states can be parametrized by the 
formula (|B"5j ) The results are presented in Fig.|5]. At small temperatures T < 80 MeV 

both the S- and P-wave give comparable contributions to the CI. At higher temperatures, the 
P-wave dominates. The D- and P-waves add small corrections to the CI at T > 140 MeV. 
The contribution of higher waves is assumed to be negligible. 

It should be mentioned that at very low temperatures, T < 30 MeV, the total CI becomes 
negative, in agreement with the results of Ref. |H| for the isotopically symmetric pion gas, 
while at high temperatures attraction dominates over repulsion. 

The exact CI is also compared in Fig.|] to various approximations widely used for the 
hadron gas analysis. It is seen that ignoring repulsion between pions overestimates the CI 
by more than 35%. On the other hand, the NRA underestimates the CI by at least 20%. 
At low temperatures T < 120 MeV both approximations become completely unreasonable. 
If one ignores both the finite resonance width and the repulsion (it corresponds to the ideal 
gas of pions and resonances) these two errors partially cancel each other. The simplest 



10 The background and higher pole terms which are present in the formulae of Ref. |^] were found 
to give a negligible contribution to the CI and are dropped in the present consideration. 
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approximation, surprisingly enough, appears to give better results than the both more com- 
plicated ones. (There remains, however, discrepancy up to about 15% at high temperatures). 
This means that both effects, the repulsion and the finite resonance width, should be taken 
into account simultaneously. Including either of these effects without another one increases 
rather than decreases the numerical errors with respect to the simplest IG model of pion 
and two-pion resonances. 

Comparing the CI 6 2 with the ideal gas CI 62 ( see Fig. |8|) one observes that the 
interactions give essential contribution to the CI already at T = 70 MeV. At T > 150 MeV 
the interaction part clearly dominates over Bose effects related to b^. 

To estimate the influence of the second CI on the thermodynamical properties of the 
pion gas we have calculated the particle density in the second cluster approximation 

n = n + 2bfz 2 , (77) 

where 

00 



n = $>J°V = -^/ dpp* — / /9 ; 9 - — (78) 



f^" 1 " 2tt 2 7o ^[Vf- 



cxp 

is the density of the ideal pion gas (without resonances). The calculations were done assum- 
ing that the chemical equilibrium, /x = 0, is reached. 

The temperature dependence of the ratio n/n is shown in Fig.|7J It can be seen that all 
approximations give consistent results up to 5 4- 15%. A rather small errors is explained by 
the fact that at low temperatures (T < 120 MeV), where ignoring either the finite resonance 
width or the repulsion between pions leads to huge errors in the value of cluster integral, the 
activity of the equilibrium pion gas is small and the contribution of the second term of the 
cluster expansion into the value of particle density is not important. On the other hand, at 
large temperatures, where the second term becomes comparable with the first one, the both 
approximation provide more exact value of the cluster integral. Again, the ideal gas model 
provides the best approximation at all temperatures, except very large ones T > 180 MeV. 

This conclussion is close to that of Ref . M , where it was pointed out that the interacting 
pion gas in the second cluster approximation only slightly differs from the ideal gas of pions 
and p-mesons due to the nearly exact cancellation of the contributions from S-wave attractive 
and repulsive channels. That is the repulsive interactions and the contribution of the broad 
cr-resonans can be dropped simultaneously. The aim of our further consideration is to 
take properly into account the hard-core repulsive interactions. In this case the a-meson 
contribution must be retained. 

It is seen from Fig.[T] that the contribution of the second term of cluster expansion to the 
particle density is comparable to that of the ideal gas. There is no reason to expect that 
the higher terms are negligible. The purpose of the next section is to go beyond the second 
cluster approximation. 

VI. VAN DER WAALS EQUATION 

Taking into account the attractive parts of higher cluster terms is straightforward: as- 
suming that the attractive interaction of three and more pions is dominated by resonance 
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interaction (as it was in the two pion case) we just add the activities of all the lightest pion 
resonances to the grand canonical partition function logarithm: 

log Z(V, T, /i) = log (V, T, /i) + V^z r , (79) 

r 

where the first term in the right hand side is the partition function of the ideal pion gas and 
the sum in the second term runs over not only the two pion resonances (Table |), but also 
includes the resonances decaying into three and more pions (Table |J). For the activity of 
each resonance species we use the following expression 

z r = r de({e)g r <t>{T;e) exp (£) . (80) 

JN r m V 1 J 

The chemical potential fi r of a resonance decaying into N r pions is proportional to the pion 
chemical potential: 

H r = N r /i (81) 

In our calculations we assume chemical equilibrium: \i r — \i — 0. For the two-pion reso- 
nances from Table | we put 

C e = --r 82 

7r de 

and use the parametrization fl64|) and (|65]), so that the contribution of two pion resonances 
is reduced to the b^z 2 , where the b^ is an attractive part of the CI b$ shown in Fig. |[ 

For the resonances from Table |lj we use the Breit-Wigner profile function fl7lD with the 
normalization (^) In the NRA the expression ( p0|) is reduced to Eq. (j62[) . 

It follows from Eq. (|79|) that the pressure and the pion density are calculated from 
the ideal gas model for pion and pion resonances with finite width (that is the repulsive 
interactions are ignored): 

P( T ^) = Po{T,n) +Tj2 Z r =P0(T,H) +J2Pr( T ^r) (83) 

r r 

n(T, n) = n (T, + N r z r = n o(T, + N rn r {T, fJ, r ), (84) 

r r 

where the ideal pion gas pressure po is given by the formula 

d 3 p , _[ ffi- VP 2 + m' 2 \] 

1 - exp ( , (85) 



Po(T^) = - 9 T J ^log 



T 



11 The formula of Ref. [Q] does not lead to large error because of the dominating contribution 
comes from the narrow resonance cj(782), for which both creteria ( |63[ ) and ( p7| ) are fulfilled. On 
the other hand, lack of detailed experimental information on the phase shifts in the vicinity on 
the broad np- and 7r<7-resonances as well as large uncertainties in their masses, width and decay 
fractions make impossible and useless the application of the more exact formula. 
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with the particle density of the ideal pion gas given by Eq. 

To take into account the hard core repulsion between pions and resonances we use the 
excluded volume Van der Waals model In the framework of this model the pressure 
pVdw ^ ^ f one-component gas of particles with the excluded- volume parameter v o can 
be found from the transcendental equation 

p VdW (T, p) = p id (T, p - v p VdW (T, p)) , (86) 

where p ld (T,p) is the pressure of corresponding ideal gas. The particle density n VdW (T, p) 
is related to that of ideal gas by the expression 



7? 



VdW 



id 

n 



(T,p-v oP VdW (T,p) 



(T ' /i) 1 + v n id (T, p - v p VdW (T, p)) ' (87) 



The above model can be straightforwardly generalized to a multi-component gas, if one 
assumes that all particle species have the same excluded-volume parameter. In our calcula- 
tions we put it to be the same for the pions and pion resonances. The standard procedure of 
derivation of the Van der Waals equation in the statistical physics shows that the excluded- 
volume parameter is equal to the absolute value of the repulsive part of the second virial 
coefficient |]TJJ. Therefore, the excluded-volume parameter for the pions can be identified 
with the repulsive part of the CI (See Fig. ||): 



(r) 



(88) 



Hence, to find the pressure of the interacting pion gas in the framework of the Van der Waals 
excluded volume model we solve the transcendental equation 

p VdW (T,p) =p (T,p) + Y,Pr(T,Pr) (89) 

r 

p = p-v p VdW (T,p) , 

Pr = Pr ~ V P [T, p) . 

The particle density of the pions is found from 

n vdw {T } = n (T p) + E r N r n r (T p r ) 

The result of the calculation is shown in Fig. |^. It is seen that essential deviation 
from the second order cluster expansion take place at the temperatures T > 140 MeV. A 
comparison with the ideal gas model of pions and pion resonances shows that the effects of 
the hard-core repulsion are not cancelled by the effects of the finite resonance width. This 
leads to an essential (up to 30%) suppression of the pion density with respect to the ideal 
gas of pions and pion resonances. 



VII. CONCLUSION 

A quantum mechanical formula for the second cluster integral for the gas of relativistic 
particles with hard-core interaction was derived and analyzed. In the nonrelativistic classical 
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limit, this formula is reduced to the expression used in Refs. []2]||. In the quantum case, 
however, the value of the cluster integral appears to be much larger in magnitude than 
the corresponding classical value and, in contrast to the classical case, depends on the 
temperature even in the nonrelativistic limit. It has been demonstrated that the second 
cluster integral for the pion gas all reasonable temperatures is far away from the classical 
limit. Its repulsive part, which can be interpreted as proper particle volume, is an order 
of magnitude larger than it could be expected from the classical evaluation. It should be 
mentioned that not only quantum effects but also relativistic ones are important in the case 
of pion gas. Surprisingly, they essentially modify the proper pion volume even at relatively 
low temperatures T ~ 30 MeV. 

The role of finite resonance width in the second cluster integral was studied. It was estab- 
lished that the widely used add hoc formula with the normalized Breit-Wigner resonance 
profile is unsuitable for broad resonances lying close to the threshold, the parametrization 
of the experimental phase shifts should be used instead. The most striking example of this 
kind is the cr-resonance. Our analysis shows that in the case of cr-resonance the calculations 
with the normalized Breit-Wigner profile can give even worse result than simple zero-width 
approximation. 

At the second order of cluster expansion, due to the presence of broad resonances in 
the 7T7r-system, the negative contribution of the hard core repulsion into the cluster inte- 
gral almost canceled by positive contributions of finite resonance widths in a rather broad 
temperature range. Because of this fact, the thermodynamical properties of the interacting 
pion gas in the second cluster approximation appear to be quite similar to those of the ideal 
gas of pions and two-pion resonances: the error in the value of the particle density does 
not exceed a few percents. Surprisingly, the account for finite resonance widths without 
account for the hard core repulsion as well as consideration of the hard core repulsion when 
the resonance widths are neglected worsen rather than improve a simple ideal gas model of 
pions and zero-width pion resonances. Both effects should be either neglected or taken into 
account simultaneously. 

This does not mean, however, that we can restrict ourselves to the simple ideal gas picture 
of pions and zero-width pion resonances at all temperatures. As it has been demonstrated, 
when the temperature is sufficiently high (T > 140 MeV), the pion density becomes so 
large that the cluster expansion cannot be truncated at the second order. In contrast to the 
second cluster approximation, an appreciable deviation from the ideal gas model is observed, 
when the higher order terms are taken into account. In the framework of Van der Waals 
excluded-volume model, the pion density appears to be up to 30% lower than that of the 
ideal pion-resonance gas. Hence, at high particle densities the correct model of the pion gas 
must include all pion resonances and the resonance width as well as the repulsive interactions 
between the particles must be taken into account. 

It should be emphasized that, if the model takes properly into account the hard core 
repulsion, there is no reason to drop the a-resonance. It must be included into the model 
along with other resonances. 

The developed in the present paper approach will be used for calculation of excluded 
volumes of other hadrons. This will allow us to study the influence of hard core repulsion 
on the properties of realistic hadron gas including (anti-)nucleons and strange particles by 
means of multicomponent Van der Waals equation M. We expect essentially larger excluded 
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volume effects for nucleons: preliminary calculations show that the proper volume of the 
nucleon is essentially (by the factor 2 -=-2.5) larger than that of the pion. Thus the hard core 
repulsion may essentially modify the particle number ratios in comparison to widely used 
ideal resonance gas model. 



ACKNOWLEDGMENTS 



We thank K.Bugaev for helpful discussions and comments. We acknowledge the financial 
support of DAAD and DFG, Germany. The research described in this publication was made 
possible in part by Award No. UP1-2119 of the U.S. Civilian Research & Development 
Foundation for the Independent States of the Former Soviet Union (CRDF). 



18 



REFERENCES 



[1] H. Stocker, W. Greiner and W. Scheid, Z. Physik A286, 121 (1978). 
H. StScker, A.A. Ogloblin and W. Greiner, Z. Phys. A303, 259 (1981). 
D. Hahn and H. Stocker, Nucl. Phys. A452 723 (1986). 
[2] G. D. Yen and M. I. Gorenstein, Phys. Rev. C59, 2788 (1999). 

P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. B465, 15 (1999) [nucl- 
th/99030iq . 



[3] G. D. Yen, M. I. Gorenstein, W. Greiner, S.N. Yang, Phys. Rev. C56, 2210 (1997). 
G. D. Yen, M. I. Gorenstein, H. Stdcker, W. Greiner, S.N. Yang, J. Phys. G24 1777 
(1998). 

[4] K.G. Denisenko and S. Mrowczynski, Phys. Rev. C35, 1932 (1987). 
[5] W. Weinhold, B. Friman and W. Norenberg, Phys. Lett. B433, 236 (1998) [nucl 
th/97100T4 . 



[6] R. Hagedorn and J. Rafelski, Phys. Lett. B97, 136 (1980). 

[7] M. I. Gorenshtein, V. K. Petrov and G. M. Zinovev, Phys. Lett. B106, 327 (1981). 

D. H. Rischke, M. I. Gorenstein, H. Stocker and W. Greiner, Z. Phys. C51, 485 (1991). 
[8] R. Venugopalan and M. Prakash, Nucl. Phys. A546, 718 (1992). 

[9] M. I. Gorenstein, A. P. Kostyuk and Y. D. Krivenko, J. Phys. G G25, L75 (1999) 
||nucl-th/ 99060681 

[10] K. A. Bugaev, M. I. Gorenstein, H. Stocker and W. Greiner, Phys. Lett. B485, 121 
(2000) |frucl-th/ 0004061 



[11] E. Beth and G.E. Uhlenbeck, Physica 4, 915 (1937). 

[12] K. Huang, "Statistical mechanics" John Wiley & Sons (1963). 

[13] B. Kahn and G.E. Uhlenbeck, Physica 5, 399 (1938). 

[14] J.E. Mayer and M. Goeppert Mayer, "Statistical mechanics" John Wiley & Sons (1977). 
[15] W. Greiner, L. Neise and H. Stocker "Thermodynamics and Statistical Mechanics" 
Springer, iVF(1995). 

[16] S. Ishida, T. Ishida, M. Ishida, K. Takamatsu and T. Tsuru, Prog. Theor. Phys. 98, 
1005 (1997). 

[17] S. Ishida, M. Ishida, H. Takahashi, T. Ishida, K. Takamatsu and T. Tsuru, particle," 

Prog. Theor. Phys. 95, 745 (1996) |hep-ph/9610325 . 
[18] S. Ishida, |hep-ph/9905260[ . 

[19] M. Abramowitz and I. Stegun, "Handbook of Mathematical Functions" Pover Publica- 
tions (1968). 

[20] B. Hyams et al., Nucl. Phys. B64, 134 (1973). 

[21] R. Dashen, S.K. Ma, H.J. Bernstein Phys. Rev. 187, 345 (1969). 

[22] A. Dobado and J. R. Pelaez, Phys. Rev. D59, 034004 (1999) [|hep-ph/9806416 . 

[23] I.N. Mishustin and W. Greiner, J. Phys. G G19, L101 (1993). 

[24] D. H. Rischke, M. I. Gorenstein, H. Stocker and W. Greiner, Z. Phys. C51, 485 (1991). 



19 



TABLES 



TABLE I. Parameters for the lightest resonances in 7T7r-system. 



Resonance 


Isospin 


Spin 


Mass 


Width 


Elasticity 


Interaction 




Ir 


Lr 


M r (MeV) 


L r (MeV) 




radius R r (GeV^ 1 ) 


a 








585 


385 






jo (980) 








993.2 


54.32 






p(770) 


1 


1 


777 


155 


1 


3.09 


/ 2 (1270) 





2 


1281 


205 


0.84 


4.94 


/O 3 (1690) 


1 


3 


1713 


228 


0.26 


6.38 



TABLE II. Parameters for the lightest resonances decaying into 3 and more pions. 



Resonance 


Degeneration 


Mass 


Width 


Elasticity 


Number of pions in 




factor g r 


M r (MeV) 


L r (MeV) 




the final state N r 


^(782) 


3 


782 


8.41 


0.888 


3 


0(1020) 


3 


1019 


4.43 


0.155 


3 


/ii(1170) 


3 


1170 


360 


~ 0.5 


3 


6i(1235) 


9 


1230 


142 


~ 1 


4 


oi(1260) 


9 


1230 


425 


~ 1 


3 


/i(1285) 


3 


1282 


24 


0.35 


4 


vr(1300) 


3 


1300 


400 


~ 1 


3 


a 2 (1320) 


15 


1318 


107 


0.70 


3 
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FIG. 1. The dependence of the CI (ro,T,m) on the temperature in the hard sphere model. 



The formula (52) is used for the relativistic calculations. The results are compared to the nonrel- 
ativistic approximation 



21 




FIG. 2. The contribution of u-resonance to (T). The exact value is compared to those 
calculated in narrow resonance approximation (NRA) and using normalized Breit-Wigner (BW) 
profile of the resonance (fT3|). 
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FIG. 3. The contribution of p(770) to b 2 (T). The exact value is compared to those calculated 
in narrow resonance approximation (NRA) and using normalized Breit-Wigner (BW) profile of the 
resonance (|73|). 



23 



4 - 



CO 







-2 - 







n 1 

- SO attractive 

SO repulsive 

- SO total 

- S2 

- S0+S2 




50 100 150 

T (MeV) 



200 



FIG. 4. The contribution of Sq and 52 states of 7T7r-scattering into the second CI b 
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FIG. 6. The exact value of the CI b 2 compared to various approximations: "No repulsion" - 
the repulsive part of the 7T7r-interaction is dropped out; "NRA" — narrow resonance approximation; 
"NRA no repulsion" — both repulsion and final resonance width are ignored. 
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FIG. 7. The ratio of the particle density of the pion gas calculated in the second order cluster 
expansion to the particle density of the ideal pion gas. Various approximations are shown: "No 
repulsion" — the repulsive part of the tttt- interaction is dropped out; "NRA" — narrow resonance 
approximation; "NRA no repulsion" — both repulsion and final resonance width are ignored. 
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FIG. 8. The attractive and repulsive part of CI and its total value. The ideal gas CI b^ is 
also shown. 
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FIG. 9. The pion density in the Van der Waals excluded volume model compared to the ideal 
gas of pions and broad pion resonances (no repulsion) , to the ideal gas of pions and zero width pion 
resonances (NRA no repulsion) and to the pion gas in the second order cluster approximation. 
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